Conservation of circulation in magnetohydrodynamics 

Jacob D. Bekenstein* and Asaf Oron^ 
The Racah Institute of Physics, Hebrew University of Jerusalem, 
Givat Ram, Jerusalem 91904, Israel 
(February 5, 2008) 

We demonstrate, both at the Newtonian and (general) relativistic levels, 
the existence of a generalization of Kelvin's circulation theorem (for pure flu- 
ids) which is applicable to perfect magnetohydrodynamics. The argument is 
based on the least action principle for magnetohydrodynamic flow. Examples 
of the new conservation law are furnished. The new theorem should be helpful 
in identifying new kinds of vortex phenomena distinct from magnetic ropes 
or fluid vortices. 

I. INTRODUCTION 

Kelvin's theorem on the conservation of circulation of a simple perfect fluid has played an 
important role in the development of hydrodynamics. For instance, it shows that potential 
flows are possible, that isolated vortices can exist, that they obey the Helmholtz laws, etc. 
Kelvin's theorem is valid only for flows in which the body force per unit mass is a gradient; 
mostly this includes incompressible or isentropic flows of one-component fluids. 

Most flows in geophysics and astrophysics are more complicated. In particular, many 
fluids in the real world carry magnetic fields: they are magnetofluids. Yet the Lorentz force 
per unit mass on a magnetofluid is almost never a perfect gradient. Thus the circulation 
theorem in its original form is almost never true in magnetohydrodynamics (MHD). Must we 
then surrender the many insights that Kelvin's theorem conferred on pure hydrodynamics ? 

Not necessarily. One might speculate that a suitable combination of fluid velocity v and 
magnetic induction B may inherit the property of having a "circulation" on a closed curve 
which is preserved as that curve is dragged with the magnetofluid. Such conserved circulation 
might play as useful a role in MHD as has Kelvin's circulation in pure fluid dynamics. For 
example, it might help characterize a set of magnetoflows as being potential in some sense, 
with consequent simplification of this intricate subject. Or it might help to characterize 
a new type of vortex, a hybrid vorticity-magnetic rope. In view of the importance of the 
vortex phenomenon in contemporary physics, this last possibility is by itself ample reason 
to delve into the subject. 

Two decades ago, E. Oron discovered, with the formalism of relativistic perfect MHD, 
a circulation theorem of the above kind. Although some of its consequences for new helicity 
conservation laws have been explored |J, this new conserved circulation has remained ob- 
scure. Contributing to this, no doubt, is the fact that it has only been derived relativistically, 
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and that this derivation is an intricate one, even for relativistic MHD. In addition, Oron's 
derivation assumes both stationary symmetry and axisymmetry, while it is well known that 
Kelvin's theorem requires neither of these. 

In the present paper we use the least action principle to give a rather straightforward 
existence proof for a generically conserved hybrid velocity-magnetic field circulation within 
the framework of perfect MHD which does not depend on spacetime symmetries. We do 
this both at the Newtonian (Sec. II) and general relativistic (Sec. Ill) levels; the importance 
of MHD effects in pulsars, active galactic nuclei and cosmology underscores that this last 
arena is not just of academic importance. 

As mentioned, we approach the whole problem not from equations of motion, but from 
the least action principle. Lagrangians for nonrelativistic pure perfect flow have been pro- 
posed by Herivel JJ, Eckart [|J, Lin 0, Seliger and Witham ||, Mittag, Stephen and Your- 
grau |7| and others. Many of the proposed Lagrangians necessarily imply irrotational flow, 
i.e. not to generic flow, a deficiency which is often missed by the authors. Lin ]5| introduced 
a device that allows vortical flows to be encompassed. This device is used by Seliger and 
Witham. Lagrangians for nonrelativistic perfect MHD flow in Eulerian coordinates have 
been proposed by Eckart ||, Henyey ||, Newcomb |I0|| , Lundgren [|ll| and others. 

In special relativity Penfield JEJ proposed a perfect fluid Lagrangian which admits vor- 
tical isentropic flow. The early general relativistic Lagrangian of Taub [H|,14] as well as the 
more recent one by Kodama et. al [15| describe only irrotational perfect fluid flows. The Lin 
device is incorporated by Schutz |16| , whose perfect fluid Lagrangian admits vortical as well 
as irrotational flows in general relativity. Carter [|l7l , |T8ll introduced Lagrangians for particle- 
like motions from which can be inferred the properties of fluid flows, including vortical ones. 
Achterberg [|n| proposed a general relativistic MHD action, which, however, describes only 
"irrotational" flows. Thompson |20j used this Lagrangian in the extreme relativistic limit. 
Heyl and Hernquist |T| modified it to include QED effects. In this paper we follow mostly 
Seliger and Witham |J and Schutz 

In Sec. II. A we propose a nonrelativistic MHD Lagrangian, and show in Sec. II. B and 
II. C that it gives rise to the correct equations of motion for the density, entropy, velocity and 
magnetic fields in Newtonian MHD. In Sec. II. D we derive from it the conserved circulation, 
defined in terms of a new vector field R, and discuss its invariance under redefinition of 
R. Sec. II. E furnishes two examples of the conserved circulation in action. In Sec. III. A 
we collect all the equations of motion of general relativistic MHD, and propose a general 
relativistic MHD Lagrangian in Sec. III.B. Sees. III.C and III.D recover all the relativistic 
MHD equations of motion from it. Finally in Sec. III.E we generalize the conserved MHD 
circulation to the general relativistic case. 



II. VARIATIONAL PRINCIPLE IN EULERIAN COORDINATES 

A. The Lagrangian Density 

Perfect MHD describes situations where the flow is nondissipative, and, in particular, 
when the magnetoflow has "infinite conductivity", and where Maxwell's displacement cur- 
rent may be neglected in Ampere's equation. We shall adopt this approximation. We work 
in eulerian coordinates: all physical quantities are functions of coordinates Xj or r which 
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describe a fixed point in space. We first summarize the MHD equations. We work in units 
for which c = 1. 

First of all, the fluid obeys the equation of continuity (d t = d/dt) 

«9 t p + V • (pv) = 0, (2.1) 

where p(r, t) is the mass density per unit volume of the fluid and v(r, t) is the fluid's velocity 
field. Second, since there is no dissipation, s, the entropy per unit mass, must be conserved 
along the flow: 

Ds = d t s + v- Vs = 0. (2.2) 

Here we have defined the convective derivative D, which in Cartesian coordinates has the 
same form for scalars or vectors. With the help of Eq. (|2.1|) this equation can be written as 



d t (ps)+V-(psv) = 0. (2.3) 

Third, "infinite conductivity" implies that E + (v/c) x B = 0, where E and B are the 
electric and magnetic fields, respectively. Combining this with Faraday's equation yields the 
so called field-freezing equation 

d t B = V x (v x B) , (2.4) 

which implies Alfven's law of conservation of the magnetic flux through a closed loop moving 
with the flow. Finally, the evolution of the velocity field is governed by the MHD Euler 
equation, 

pJv^-Vp- P VC/ + < VxB > XB , (2.5) 

where p is the fluid's pressure (here assumed isotropic), and U(r,t) is the gravitational 
potential. 

The least action principle is in general 

SS[f a ] =5jdtjd 3 r C(f a , d t f a , V/ a ) = 0. (2.6) 

Here the action S is a functional of various fields f a (r,t), a = 1,2, ■ ■ -. One varies each 
f a , transfers time and space derivatives of each variation 5f a to the adjacent factor by 
integration by parts, and sets to zero the overall coefficient of the bare 5f a . This gives us 
the Lagrange-Euler equation 

It is usually more convenient to get the equation for each f a ab initio by the above procedure, 
rather than by using Eq. ( |2.7|) . 

We now propose the following Lagrangian density for MHD flow of perfect infinitely 
conducting fluid which incorporates Eqs.( |2.1| - |2.4j ), as three Lagrange constraints 
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£ = pv 2 /2 - pe (p, s)- P U - B 2 /(8tt) + 
+ <j) [d t p + V ■ \pv)] + rj [d t (ps) + V ■ (psv)] 

+ A [dt(jry) + V • (ptv)] + K • [9 t B — V x (v x B)] . (2.8) 

In the above e (p, s) is the thermodynamic internal energy per unit mass; in the total La- 
grangian the corresponding total internal energy enters as a potential energy. The magnetic 
energy, the volume integral of B 2 / (8tt), also enters the total Lagrangian as a potential energy. 

In Eq. ( |2.8|) 0, i] are Lagrange multiplier fields which locally enforce the conservation 
laws ( |2.1| - |2.2p , as may be verified by varying with respect to these multipliers. K is a triplet 
of Lagrange multiplier fields which enforce the field-freezing constraint Eq. ( |2.4j ): varying 
with respect to K reproduces Eq. (|2.4j) at every point and time. Finally, A is a Lagrange 



multiplier field which enforces the Lin constraint on a new field, 7: 

dt(jry) + V • (pv 7 ) = or £> 7 = 0. (2.9) 



Here we have used Eq. (|2.1|) to reduce to the second form. Lin's field 7, like s, is conserved 
along the flow, but unlike s it does not occur elsewhere in the Lagrangian. Lin interprets 
7(r, t) as one of the three initial Lagrangian coordinates which label each fluid element. But 
whatever the interpretation, the condition ( |2.9|) is essential so that the flow can be vortical 
also in the limit B — > 0. This matter is further discussed in the following section. 



B. The Equations of Motion 

Can our proposed Lagrangian density reproduce all the equations of motion of perfect 
MHD flow ? We have already seen that it does reproduce Eqs. ( f2~l~] - [2.2D and ( |2.4|) . Let us 

now vary 7 to get 

D\ = 0, (2.10) 

so that A, like 7, is conserved with the flow. Both this and Eq. (|2.9|) will be essential in 
demonstrating the existence of the new conserved circulation. Next we vary s; remembering 
that (de/ds) p is just the fluid's temperature T, we have 

D V = -T, (2.11) 

which establishes that 77 decreases along the flow. The next variation is one with respect to 
p. Recalling that (de/dp) s = p/p 2 , introducing the enthalpy per unit mass w = e + p/ p, and 
using Eqs. (|2 . 1 0| - |2 . 1 1]) we get 

D<p = v 2 /2- w + T-U. (2.12) 

When we vary v in the action we may take advantage of the identity V • (A x B) = 
B-VxA — A-VxB and Gauss' theorem to flip the curl operation from 5v x B onto K. 
Then the identity A ■ B x C = — B • A x C helps to shift the <5v into the position of a factor 
in a scalar product. We may then factor out the common 5v and isolate the vector equation 

v = V0 + 7VA + sVr? + Q (2.13) 
QEBxR/p. (2.14) 
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where R = V x K. This is neither a solution for v (A and rj not known), nor an equation of 
motion (v appears undifferentiated). In the next subsection we show that this prescription 
for v leads to the MHD Euler equation ( |2.5| ). 

Expression Q2.13Q shows the importance of including Lin's field 7. For suppose we con- 
sider an unmagnetized fluid in isentropic (s = const.) flow. Without 7 the expression for 
v is a perfect gradient, which means the proposed Lagrangian density describes only irrota- 
tional flows, a small subset of all possible ones. It is well known |],[7| that this problem does 
not appear when one couches the problem in Lagrangian coordinates because one gets then 
an equation, not for v, but for the fluid's acceleration. Lin's ]5[] way out of this difficulty 
is to remember that the initial coordinates of the fluid element are maintained throughout 
its flow. These coordinates "label" the element, and this can be interpreted as a triplet of 
constraints (one for each coordinate) of the form Aj (dbi/dt + V&j) , where b is the initial 
vector coordinate for the element in question. Lundgren used this triplet form for the 
MHD case. It was later shown (see for example |J) that the triplet can be reduced to a 
single constraint with the help of Pfaff 's theorem. One thus returns to form ( |2.8| ) of the 
Lagrangian density and Eq. ( |2.13| ) for the fluid velocity. The vorticity is now (still excluding 
B) 

w = Vxv = V7xVA + VsxV)j, (2.15) 

so we see that isentropic vortical flow is possible. 

In the MHD case, the magnetic term in Eq. ( |2.13| ) contributes to the vorticity. Henyey 
||, who suggested a Lagrangian density similar to ours, occasionally dropped the Lin term 
in the MHD case. However, we shall retain the Lin term throughout. It might seem peculiar 
at first that adding a constraint like Lin's permits the appearance of solutions (vortical) 
which were forbidden before it was imposed. But we must remember that we add to the 
Lagrangian not only a constraint, but also a new degree of freedom, j(r, t), and it is natural 
that with more degrees of freedom the class of allowed flows will expand. 

Finally, we vary B in the action; by similar manipulation to those which gave Eq. Q2.13p 
we get 

d t K = v x R- B/(4tt). (2.16) 

Taking the curl of this equation we get the more convenient one 

d t R = V x [v x R - B/(4tt)] = V x (v x R) - J. (2.17) 

Here J = V x B/47r is the electric current density coming from Ampere's equation. Notice 
the similarity between Eq. ( J2.17 ) and ( |2.4| ). Eq. ( |2.17[ ) says that the rate of change of the 



flux of R through the surface spanning a closed curve carried with the flow equals minus 
the flux of the electric current density through that curve. 



C. The MHD Euler equation 

We now show that the Lagrangian density ( |2.8| ) yields the correct MHD Euler equation. 
We first operate with the convective derivative D on Eq. (|2.13|) remembering that Ds = 
and -D7 = 0: 
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Dv = DV(p + 7DVA + sDVi] + DQ. (2.18) 
We now use the identity 

DV = WD — (Vv) • V, (2.19) 

where in Cartesian coordinates 

[(vv).v],,x;gA, ( , 20) 

j J 

in conjunction with Eqs. (|2 . 1 0| ~ pTT2f) to transform Eq. (|2.18|) into 

Dv = V{v 2 /2 -w + Ts-U)- sVT - s(Vv) • V77 

- (Vv) • V</> - 7(Vv) ■ VA + DQ. (2.21) 

From the thermodynamic identity dw = Tds + dp/ p we infer 

-Vw + TVs = -Vp/p, (2.22) 

and we also have Vv 2 /2 = (Vv) • v, where the meaning of the right hand side is clear by 
analogy with Eq. ( p720|) . Thus Eq. ( |272T| ) turns into 

Dv = -Vp/p -VU + (Vv) • (v - V0 - sVr] - 7VA) + DQ. (2.23) 

Finally comparing with Eq. ( f2.13j ) we see that the last brackets stand for Q so that 

= -Vp/p -VU + (Vv) ■ Q + DQ. (2.24) 



Thus, magnetic term aside, we have recovered the Euler equation (|2.5| ). 

We now go on to calculate the Q dependent terms. We may rewrite the equation of 
continuity ( [2.1|) as 



Dp = -pV • v, (2.25) 

With this, the Gauss law V ■ B = and the identity V x (A x B) = B • VA - BV ■ A - 
A • VB + AV • B, Eq. (|2.4|) may be recast in the well known form 



D(B/p) = ((B/p) ■ V) v. (2.26) 
Analogously, because V • R = 0, Eq. ( |2.17| ) may be put in the form 

D~R = (R • V)v — RV v — J. (2.27) 

Therefore, 

(Vv) • Q + DQ = -(Vv) • (R x B/p) — DR x B/p — R x D(B/p) 
= -(Vv) • (R x B/p) - ((R • V)v) x (B/p) 

+ (V • v)R x (B/p) - R x ((B/p) • V) v + J x B/p. (2.28) 

The four terms in the second version of Eq. (|2.28|) involving derivatives of v can be 
shown to cancel out by expanding them out in Cartesian coordinates. Hence, Eq. ( |2.24| ) is 
the magnetic Euler equation with the usual Lorentz force per unit mass, J x B/p, in addition 
to the pure fluid terms. The fact that we obtain the correct MHD equations ( [2.1| - [2.2| ), ( |2.4| ) 
and ( p.5|) is testament to the correctness of our proposed Lagrangian density Eq. (|2.8|) . Note 
that Lin's field 7 has disappeared from the final equation of motion. 
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D. Circulation Conservation Law 



With the help of the above formalism, we can now prove the existence of a generalization 
of Kelvin's circulation theorem applicable to perfect MHD. Let us calculate the line integral 
of the vector 

Z = v + RxB/p (2.29) 
along a closed curve C drifting with the fluid: 

T = iz-dr. (2.30) 



ic 

According to Eq. Q2.13|) this integral is 



r = I V0 ■ dr + I 7V A -dr+i sVt? ■ dr. (2.31) 

J C J C J c 

The term involving obviously vanishes (we assume all the Lagrange multipliers are single 
valued). For like reason so does the term involving rj in the isentropic (s = const ) case 
as s can be taken out of the integral. The middle integral can be written § c 7 dX, where 
dX = VA ■ dr. But Eqs. ( |2.U| - |2.1(J| ) tell us that both 7 and A are conserved along the flow. 
Hence T remains constant as C drifts along with the flow. Since in the limit B — > 0, T 
becomes Kelvin's circulation, we have found an extension of Kelvin's theorem to perfect 
MHD. Obviously the conservation of T implies the conservation of the flux of V x Z through 
C. 

The vector field R is not unique for a given physical situation. For example, the change 
R — > R + /cB (k a real constant) leaves invariant all equations of motion, Eqs. ( |2.9| - p. 14]) , 



( 2.17 ), and ( |2.24| ), as well as the conserved circulation expressions ( |2 . 29j - pT30|) . In addition 



suppose that at time t = we define an arbitrary solenoidal (divergence-free) field b all 
over the flow, and then evolve it in time as a passive vector, i.e., in accordance with the 
frozen-in field equation (2li). Comparing with Eq. ( 2.1 7| ) we see that R+A;b and R obey the 



same equation, and both are permanently solenoidal [this property is obviously preserved 
by Eqs. ( |2.4|) and ( 2.17 ) in the MHD approximation]. 



If in Z we use R + kh in lieu of R to construct the conserved circulation, T gets the 
additional contribution 

Ar = k <£ (b x B/p) • dr = k I B • (dr x b/p) (2.32) 
Jc Jc 



Here we have used a well known vector identity. Now by analogy with B, b obeys Eq. (p. 26 ) 



which tells us that any two elements of the fluid permanently lie on one and the same line 
of b/p, and their distance, if small, is proportional to |b|/p P2|| . We can always make b 



small. Then dr x b/p is a vectorial element of area of a narrow closed strip carried along 
by the fluid, one of whose edges coincides with C. The integral in Eq. ( |2.32| ) is just the flux 



of magnetic induction through this strip (not through the space bounded by the strip), and 
we know this is conserved by virtue of Alfven's law. 

Thus with the change R — > R + kb we added some conserved magnetic flux to T, and 
did not get a new conserved circulation. The MHD flow {B, v,p,p} is evidently unchanged 
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because the MHD Euler equation (|2.5|) does not contain R, so we must conclude that in 
the expression for v, Eq. ( 2.13 - 2.14 ), the change of the Q term must be compensated by 
suitable changes in the Lagrange multipliers + sr] and A (recall that we are working with 
s = const.). Indeed, the initial choice of b involves a choice of two functions because of the 
V • b = constraint, so that the two functions <fi + srj and A are just enough to absorb the 
change R — > R + kb thus generated and leave v unchanged. It is not possible to eliminate 
R altogether by the change R — > R + kb because R and b obey different equations. This 
means the circulation conservation law we have found cannot be reduced to an Alfven type 
law; it is a new law. 

In Sec. [lllEj we shall discuss the freedom inherent in R by a covariant procedure. Fixing 
the freedom is a necessary step in any attempt to exhibit explicitly the conserved circulation. 



E. Examples 

First consider a situation where the fluid is isentropic but not flowing: v = 0. It follows 
from Eq. ( [2.1|) that p = po(r), and from Eq. (|2.4j) that B = B (r). From these facts and 
Eq. fl2.17|) we see that 

R= -t V x B (r)/(4vr) + R (r). (2.33) 

Although the physical quantities are stationary, R is not. This is so because like the elec- 
tromagnetic potential, R is not a measurable quantity, being subject to "gauge changes" 
R — > R + b as already discussed. According to Eq. fl2.29| ) the conserved circulation (around 
a contour fixed in space because v = 0) should be 

r = _ ,(VxB )xB dt+ ,R^B, ^ 
Jc 4-7rp Jc 47rp 

On the face of it, the time dependence of the first term in this simple situation puts the 
claimed circulation conservation law in jeopardy. However, according to the magnetic Euler 
equation (|2.5|), the first integrand here is equal to VU + Vp/po which, by virtue of Eq. ([2.1 9|) 
and the isentropic nature of the fluid, is a perfect gradient (for isentropic flow Vp/po = Vw). 
Hence the first integral vanishes, and the circulation is indeed time independent as required 
by our theorem. 

As a second example consider an axisymmetric differentially rotating fluid exhibiting a 
purely poloidal magnetic field. Let the flow also be isentropic and stationary. We choose 
to work in cylindrical coordinates {g, <f), z}; the hat symbol will denote a unit vector in the 
stated direction. It then follows that p = po(r), B = B (r) and v = Qg<p, where fl(g,z) 



is the angular velocity of the fluid. It is well known fl2"3| , |24|| that for axisymmetric fields 
the curl of a poloidal field is a toroidal one, and the toroidal field has only a component. 
Therefore, the electric current density J = V x B/(47r), is everywhere collinear with v 



and time independent. Since the problem is stationary, f2 satisfies Ferraro's [25,24| law of 
isorotation B • Vf2 = 0. In addition the field must be torque-free ]24| i.e. no Lorentz force in 
the (f) direction. This condition is identically satisfied for a purely poloidal field. Combining 
all of the above we get the following solution of Eq. ( |2.17|) : 

R = -tJ (2.35) 



s 



According to Eq. Q2.29Q the conserved circulation should be 



j>VLg 2 



fJxB , , 

t i — • dr 2.36 

Jc 4vrp 



where we have exploited the axisymmetry to rewrite the first term. We now verify that this 
circulation is indeed conserved. Because of the differential rotation, the contour C is grad- 
ually deformed in the azimuthal direction. The difference dtp in the azimuthal coordinates 
between two infinitesimally close fluid elements lying on C can be written as dip = d<po + 1 dQ 
where d<po is the initial difference in azimuthal coordinates while dfl is the difference between 
the elements' angular velocities. Hence we have 

I Qg 2 d(p = I Qg 2 d(p + t I Qg 2 dQ. (2.37) 

•J c J C J c 

Note that the first term is time independent while the second one is linear in time. 
The magnetic Euler equation ( |2.5| ) in cylindrical coordinate reads 



, 2 s Vp „ rr , J x B 



-il 2 gg = ^-W + -. (2.38) 

Po 4vrpo 

Again, by the isentropic condition we can write Vp/po = Vu>. Taking the integral round C 
of both sides of Eq. (|2.38|) we have 



/ n 2 gdg = I ■ dr. (2.39) 

JC JC 4:7100 



k 47rpo 

Substituting from Eq. ( ggg ) and Eq. ( g37p into Eq. (g3J ) we get 

r = I ttg 2 d(p + t I ttg 2 dtt + t <L Vt 2 gdg 
= ^ Qg 2 d<p + \j c d(p?g 2 ) = £ Qg 2 d<p , (2.40) 

and T is indeed time independent. Note that it is possible to add to R in Eq. (|2.35|) 
an arbitrary time independent solenoidal vector field Ro(r) which satisfies Ro x v = Vx- 
However, as already stressed in the previous subsection, this will only add to T a time 
independent quantity. 

It is important to note that although the example specifically relates to an axisymmetric 
problem, Eq. ( [2.35|) applies to all stationary MHD flows which have J collinear with v. 
Accordingly, T will be conserved in all such flows. 



III. RELATIVISTIC VARIATIONAL PRINCIPLE 

In this section we formulate a Lagrangian density for MHD flow in the framework of 
general relativity (GR). Greek indices run from to 3. The coordinates are denoted x a = 
(x°, x 1 , x 2 , x 3 ); x° stands for time. A comma denotes the usual partial derivative; a semicolon 
covariant differentiation. Our signature is {—,+,+, +}. We continue to take c = 1. 
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A. Relativistic MHD Equations 



The general relativistic (GR) equations for MHD were formulated by Lichnerowicz [EB 



Novikov and Thorne [27[], Carter |T7|], Bekenstein and Oron |L| and others. The role of the 
mass conservation equation (|2.1| ) is taken over by the law of particle number conservation, 

N a ;a = {nu a ). a = 0, (3.1) 

where N a is the particle number 4-current density, n the particle proper number density 
and u a the fluid 4-velocity field normalized by u a u a = —1. If s represents the entropy 
per particle (not per unit mass as in Sec. II), then the assumption of ideal adiabatic flow, 
Eq. (|2.2p, can be put in the form 

(sN a ). a = or u a s a = 0. (3.2) 

Because the flow is assumed adiabatic, the energy momentum tensor for the magnetized 
fluid is that of an ideal fluid augmented by the electromagnetic energy-momentum tensor: 

T <*0 = pg *P + ( p + p ) „V + (F a ~<F% - -F^F yS # a/3 )/(47r). (3.3) 

4 

Here p represents the fluid's energy proper density (including rest masses) and p the scalar 
pressure (again assumed isotropic), while F a/3 denotes the electromagnetic field tensor. As 
usual the covariant divergence T a/3 ;j g must vanish (energy-momentum conservation). In 
consequence T al3 .p + u a u 1 T 1 ^ \p = which can be recast as 

(p + P Vu a . fi = -(g°0 + u a u^ + F^F^.^/^71). (3.4) 

The term a a = u^u a .^ stands for the fluid's acceleration 4-vector. The effects of gravitation 
are automatically included by the appeal to curved metric and covariant derivatives. This 
equation parallels Eq. (|2.5|) ;as usual in GR the pressure contributes alongside the energy 
density to the inertia. The electromagnetic field tensor is subject to Maxwell's equations 

F^.n = 4vrT (3.5) 

F<xP,>y + Fpy,a + F-ya,[3 = 0. (3.6) 

where j a denotes the electric 4-current density. Putting all this together we have the GR 
MHD Euler equation 

( p + p)a Q = -^% + F^ J/3 , (3.7) 
where we have introduced the projection tensor 

The above equations do not completely specify MHD flow (as opposed to flow of a generic 
magnetofluid) . For any flow carrying an electromagnetic field, the (antisymmetric) Faraday 
tensor F a p may be split into electric and magnetic vectors with respect to the flow: 
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E a = F af) vP (3.9) 
B a = *F /3a u (3 = l -e Pcn5 i^V. (3.10) 

Here e a ^s is the Levi-Civita totally antisymmetric tensor (eoi23 = (— g) 1 ^ 2 with g denoting 
the determinant of the metric g a @) and *F a! 3 is the dual of F a p. In the frame moving with 
the fluid, these 4-vectors have only spatial parts which correspond to the usual E and B, 
respectively. The inversion of Eqs. Q3.9H3.10Q is 

F af3 = u a Efi - u p E a + e^su'B 5 (3.11) 

For an infinitely conducting (perfect MHD) fluid, the electric field in the fluid's frame must 
vanish, i.e., 

E a = F aP u p = 0. (3.12) 
This corresponds to the usual MHD condition E + v x B = 0. 



B. Relativistic Lagrangian Density 



Inspired by Schutz's [|nj Lagrangian density for pure fluids in GR, we now propose a 
Lagrangian density for GR MHD flow which reproduces Eqs. ( |3.1| - |3.2| ), ( p. 5^3.7] ) and ( |3.12| ). 



Like Schutz we include Lin's term, which proves essential to our subsequent proof of the 
existence of a circulation theorem. The proposed Lagrangian density is 

C = -p(n, s) - F a pF a ?/(lQn) + 0iV% + V (sN a ). a + A ( 7 iV Q ). Q + r a F a0 N^. (3.13) 

Now in GR the scalar density C (— g) 1 / 2 replaces C in the action (|2.6|), and is what enters in 
the Euler-Lagrange equations ( |2.7|) . The covariant derivatives cause no problem; for example 
{—g) 1 ^ 2 (j)N a . a = (f>[(— g) l l 2 N a ] y0l) whose variation with respect to N a is easily integrated by 
parts. 

As in the nonrelativistic case, <p is the Lagrange multiplier associated with the con- 
servation of particle number constraint, Eq. ( |3.1| ), rj is that multiplier associated with the 
adiabatic flow constraint, Eq. (|3.2|) , and A is that associated with the conservation along the 
flow of Lin's quantity 7. We view 7, N a and s as the independent fluid variables, while n 
and u a are determined by the obvious relations 

- N a N a = n 2 ; u a = n' 1 N a . (3.14) 

Strictly speaking, one should include in £ a new Lagrange multiplier times the constrained 
expression N a N a + n 2 . Rather than clutter up £ further, we enforce this constraint below 
by hand. 

As usual, we view the vector potential A a , rather than the electromagnetic field tensor 
F a /3 = Afo a — A a .p = Ap >a — A a> /3, as the independent electromagnetic variable. In con- 
sequence, the Maxwell Eqs. ( |3.6| ) are satisfied as identities. The last term in C enforces 
the "vanishing of electric field" constraint, Eq. ( |3.12| ); r Q is a Lagrange multiplier 4-vector 



field. Because here we enforce the "vanishing of electric field" rather than the equivalent 
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flux freezing condition ( |2.4| ), the r a is more like R of Sec. II. B. than like K. Not all of r a is 
physically meaningful. For suppose we add an arbitrary function f(x") multiplied by N a to 
T a . This increments the Lagrangian density by fN a F a pN^ which vanishes identically by the 
antisymmetry of F a p. So r a and r a + f'N a are physically equivalent. We shall exploit this 
to substract from r a its component along u a . So henceforth we may take it that r a u a = 0. 

Much freedom is still left in r a . Suppose we add to it a term proportional to n~ l B a . By 
Eqs. ( |3.9[ - |3.11| ), this adds to the Lagrangian density the term E a B a . Of course we cannot 
take this to vanish at the Lagrangian's level because we have not yet obtained the freezing-in 
condition (3.12) from it. However, it is known that B a E a = ^e^^F^F^s- By introducing 

the potential A a we can write this as \ e a ^ y5 F a gA~ — ^^^ s F a B-sA~, where we have used 



:S 



the fact that e 01 ^ 5 has vanishing covariant derivatives. Obviously the last term vanishes 
by virtue of Maxwell's equations Q3.6|) which are identities in the present approach. When 
multiplied by (— g) 1 ^ 2 , the first term becomes a perfect derivative. Such term, when added 
to the integral forming the Lagrangian, is known not to affect its physical content. Thus 
r a and r a + const, x rT l B a are physically equivalent, and this transformation respects the 
condition r a u a = because B a u a = [see ( 3.10| )1. However, there is not enough freedom 
in the constant to allow us to eliminate the component of r a along B a . But in Sec. [III E| we 
shall exploit what we have just found. 



C. Equations of Motion 

We can now derive the equations of motion. Variation of <fi recovers the conservation of 
particles N a ;a . Variation of A with subsequent use of the previous result yields 



0. 



If we vary 7 we get 



X n u a = 0. 



(3.15) 



(3.16) 



These two results are precise analogs of Eqs. fl2.9|) and (12.10Q ; they inform us that 7 and A are 
both locally conserved with the flow. In view of the thermodynamic relation n~ 1 (dp/ds) n = 
T, with T the locally measured fluid temperature, variation of s gives 



u 77 



—T\ 



(3.17) 



this is the analogue of Eq. ( |2.11|) . 

We now vary N a using the obvious consequence of Eq. ( |3.14| ) , 

5n = -u a 5N a , 



together with the thermodynamic relation |27| involving the specific enthalpy 

fi = (dp/dn) s = n' 1 (p + p); 
we get the GR analog of Eq. ( |2.13| ) , 



(3.18) 



(3.19) 
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(3.20) 



The importance of Lin's 7 is again clear here; in the pure isentropic fluid case (_F Q/3 = and 
s = const.), the Khalatnikov vorticity tensor given by 



u aP = (jiup) - (fiu a ) = (7 A,/?) - (7\ a ) 



3 



(3.21) 



would vanish in the absence of 7, thus constraining us to discuss only irrotational flow. 

By contracting Eq. ( |3.20D with u a and using u a u a = — 1 as well as Eqs. Q3.12Q and 
( P^3T7|) , we get the following GR version of Eq. fl2~T2l) : 



Ts. 



(3.22) 



Thus the proper time rate of change of <fi along the flow is just minus the specific Gibbs 
energy or chemical potential. The apparent difference between the result here and Eq. (|2.12| ) 
stems from the fact that proper time rate (here) and coordinate time rate (there) differ by 
gravitational redshift and time dilation effects. These effects are not noticeable when one 
compares Eqs. (|3.17f ) with (|2.11| ) because the first refers to locally measured temperature 
and the second to global temperature; these two temperatures differ by the same factors as 
do proper and coordinate time. 

Turn now to the variation of A a . Because of the antisymmetry of F a p, the last term 
of the Lagrangian, Eq. ( |3.13|) , can be written as {r l3 N a — T a N^)A a ^. The variation of 
A a in the corresponding term in the action produces, after integration by parts, the term 
[(—g) 1 / 2 {r a N 13 — T (3 N a )] pbA a . Because for any antisymmetric tensor t a/3 , (— g) 1 / 2 ^^ = 
[(— gY^ 2 t a/3 }^, this finally leads to the equation 



4tt (V 



(3.23) 



Comparison with Eq. (3J5) shows that this result gives us a representation of the electric 
current density j a as the divergence of the bivector T a N l3 — r^N a . Such representation makes 
the conservation of charge (j a - a = 0) an identity because the divergence of the divergence of 
any antisymmetric tensor vanishes. This equation is the GR analogue of Eq (|2.17| ). Formally 
Eq. ( |3.23| ) determines the Lagrange multiplier 4- vector r a , modulo the freedom inherent in 
it. 



D. MHD Euler Equation in General Relativity 

Our central task now is to show that the equations in Sec. III.C lead uniquely to the 
GR MHD Euler equation (|3.7|) . We begin by writing the Khalatnikov vorticity uop a in two 
forms, 

U(i a = fl,pU a - /I }a Up + /lUatf - /lUp-a, (3.24) 

as well as by means of Eq. ( |3.20| ) 

+ T S , p F aS - T S , a F p6 + r 5 F a5;/3 - T S F ps , a . (3.25) 
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Contracting the left hand side of the first with N a , recalling Eq. ( |3.14j ) and that by normal- 
ization u a u a .f3 = whereas vPu^ = a a , the fluid's 4-acceleration, we get 

up a N a = —n\ifi — np^ a u a up — npap = —nhp a p^ a — npap. (3.26) 



Now contracting Eq. ( |3.25| ) with N a and using Eqs. ( |3.15| - |3.l7| ) and ( |3.12| ) to drop a 
number of terms we get 

u Pa N a = -nTs, p - T 5 , a Fp S N a + r 5 F a&]fi N a - r 5 F p5 , a N a . (3.27) 



By virtue of Eq. (p.2|), —nTs^ is the same as —nThp a s >a . It is convenient to use the 
thermodynamic identity dp = n _1 dp + Tds, which follows from Eq. ( |3.19| ) and the first law 
d(p/n) = Tds — pd(l/n), to replace in Eq. (|3.27|) —nTs^ by hp a {—np a + p >a ). Equating 
our two expressions for ujp a N a gives, after a cancellation, 



- (nfiap + hfp^) = -T d , a F ps N a + r d F aS , p N a - T d F pS]a N a . (3.28) 

The last two terms in this equation can be combined into a single one by virtue of 
Eq. ( |3.6| ), which, as well known, can be written with covariant as well as ordinary derivatives. 
Further, by Eq. ( |3.19| ) we may replace rip, by p + p. In this manner we get 



(p + p) a p = -h p a p, a + F Pa , S T 5 N a + F p& T 5 , a N a . (3.29) 

The term r s . a N a here can be replaced by two others with help of Eq. (|3.23|) if we take into 
account that N^-^ = 0: 

(p + p)a p = -hp a p, a + F ps F Sa . a /(47r) + F Pa]S T 5 N a + F ps (r a N s ) ^ . (3.30) 

We note that the last two terms on the right hand side combine into {^Fp a N a T S ^j which 

vanishes by Eq. ( |3.12| ). Now substituting from the Maxwell equations ( p.5| ) we arrive at the 
final equation 

(p + p) a/3 = -hp a p )a + Fpsj 5 , (3.31) 



which is the correct GR MHD Euler equation (|3.7|) . We have not used any information 
about r a beyond Eq. Q3.23Q ; hence Euler's equation is valid for all choices of r a . Since we 
are able to obtain all equations of motion for GR MHD from our Lagrangian density, we 
may regard it as correct, and go on to look at some consequences. 

E. New Circulation Conservation Law 



Eqs. (|3.20|) and ( |3.15| - |3.16l ) allow us to generalize the conserved circulation of Sec. II. D 



to relativistic perfect MHD. Let T be the line integral 



T = j c z a dx a , (3.32) 



where C is, again, a closed curve drifting with the fluid, and 
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Z a = jlU a — T^Fap. (3.33) 

According to Eq. ( |3.20| ), z a = a + s?7 )Q + 7A )Q . Since iQ is a gradient, its contribution 
to T vanishes. Likewise, for isentropic flow (s = const.) the term involving si]^ a makes no 
contribution to T. Thus 

r -f\ Q dx a = i jdX. (3.34) 

c ' Jc 



By Eqs. ( |3.15| - |3!T6| ) both 7 and A are conserved with the flow. Thus T is conserved along 



the flow. Note that by virtue of 7's presence, T need not vanish. 

In the absence of electromagnetic fields and in the nonrelativistic limit (/i — > m where 
m is a fluid particle's rest mass), T for a curve C taken at constant time reduces to Kelvin's 
circulation. On this ground our result can be considered a generalization of Kelvin's circula- 
tion theorem to general relativistic MHD. We have gone here beyond Oron's original result 
U in that no symmetry is necessary for the circulation to be conserved. 

To manifestly exhibit the conserved circulation, one has to know r a explicitly. The first 
step is to understand the freedom left in r a beyond that discussed in Sec. [Hi B| . The second 
is to determine r a in a specific flow exploiting for this the symmetries and other information. 
Below we address the first step; the second is left mainly to future publications. 

Given a specific MHD flow as background, let us at define a generic test field f a p = —fp a 
which satisfies Maxwell's homogeneous equations ( |3.6|) as well as the freezing-in condi- 
tion (|3.12|) , e.g. e a = fapu 13 = 0. We think of f a p as very weak, so that it does 



not disturb the MHD flow or the spacetime geometry; it is a passive tensor. Because 
fapu 13 = 0, fa/3 has only three independent components. Therefore, its full content is 
reflected in the "magnetic 4-vector" b a = ^e^a-ys f^u^, which is obviously orthogonal 
to u a . The transformation r a — > r a + kn~ l b a (k a real constant) is not a symmetry 
of the Lagrangian. However, it does not disturb the inhomogeneous Maxwell equations 
Q3.5| . |3.23|) . This is because the change in r a merely adds to the electric current the term 



b a u^ — b^u ) = (—g) x l 2 {—gY^ib 01 ^ — bPu a ) . Because of the condition e a = 0, we 
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1/2 aPyS 



may easily invert the analog of ( |3.11|) to get b a vP — b^u a = \e a,3 ' yS f^s- But since (— gY^t 
is just the constant antisymmetric symbol, our assumed equations f a p n + f ia ,/3 + fi3-y,a = 



imply that (b^u 13 — b^u a j = so that the Maxwell equations ( |3.23|) are invariant under 

r Q — *• r a + kn' 1 b a . So is the Euler equation since its derivation used only the information 
about r a inherent in ( |3.23|) . 

The expression for u a , ( |3.20|) , does seem to change under r Q — > r a + krT 1 b a , and we also 
note that r -> T + k § n~ l b p F a0 dx a . Now since the "magnetic 4-vector" b a is frozen in, like 
any such infinitesimal field, it evolves in such a way that rT x b a gives for all time that part 
of the spacetime separation of two neighboring fluid elements which is orthogonal to u a [|TJ , 
c.f. discussion after Eq. (|2.32|) . Thus rT l b a can be employed to define a thin closed strip 



dragged with the fluid such that one of its edges coincides with the curve C. Therefore, the 
increment § n~ l b^F a p dx a is just the conserved magnetic flux through this strip. Evidently 
the transformation r a — > r a + kn~ x b a has not changed the nature of the conservation law 
for T, but only added a trivially conserved quantity to it.. 

Now the MHD flow {B a ,u a ,n, p, fi} is evidently unchanged because neither the MHD 
Euler equation (p.4|) nor Maxwell's equations were changed, so we must conclude that in 



15 



the expression for u a , Eq. ( |3.20|) , the change of the r^F a p term must be compensated by 



suitable changes in the pair of Lagrange multipliers (j) + sr] and A (since we are assuming 
s = const.). They are capable of this because b a has only two independent components. 
For the condition b a u a = eliminates one of the four. In addition b a comes from f a p which 
satisfies Eqs. (|3.6|); in particular, /i 2j 3 + ^31,2 + /23,i = in the chosen coordinates. But 
since no time derivatives appear in it, this last equation serves as an initial constraint on 
b a just as the Gauss equation V • B = does for the true magnetic field. Accordingly, one 
further relation exists between components of b a so that the generic b a contains only two 
free functions. Thus the change in r^F a p can be taken up by changes in the two functions 
4> + si] and A so that fiu a is unchanged. 

Note that it is not possible to "get rid" of r a by means of the transformation r a — > 
r a + kn~ l b a because, as we shall make clear presently, r a and b a obey different equations 
of motion. Thus there must be a residual part of r a which is not changed by the trans- 
formations. It is this part which is responsible for the conserved circulation, so that the 
conservation of T cannot be reduced to magnetic flux conservation. 

The following algorithm can be used to find r a . Maxwell's inhomogeneous equations 
( p.23| ) which say that the divergence of a certain tensor vanishes can always be "solved" by 
the prescription 

F°0 - Aix{r a N p - r p N a ) = X -e a ^ 6 T 1& (3.35) 

where the new field Tys just has to satisfy Maxwell's homogeneous equations ( p.6|) , i.e. 
Tys = As,j — Ay,5. Taking the dual of Eq. ( |3.35| ) with help of the identity e-ySaP e = 
-2(5/5/'- 5J5f) gives 

*F j5 - Ane lSaP T a NP = -T 1& (3.36) 

Contracting this equation with u 1 gives the further requirement on T a p. 

T h u< = B s , (3.37) 

where we have used Eq. ( |3.10| ). The J-'s-y can always be solved for: because of gauge freedom 
there are three independent components in A a , and this is enough to find a solution for an 
arbitrary field Bs obeying B a u a = (thus three components at most). If fact, Bs does not 
determine J-'s-y uniquely: if one adds to this last one of the frozen-in f^s we discussed earlier 
in this section (which also satisfy the homogeneous Maxwell equations), Eq. (|3.37|) is still 
satisfied because fsyU 1 = 0. 

We get r a by contracting Eq. (|3.35| ) by up and remembering that F al3 U/3 = and r^up = 
0. Thus 

r a = (%Tm)- l e a ^ l5 F 1 sUp (3.38) 

It is interesting that Bs plays the role of electric part of Tsy while r a enters like the magnetic 
part of this tensor, c.f. Eq. (|3.10|) (but because Ts^u 1 7^ 0, r a evolves differently from a 
magnetic type field like B a or the b a ). It should also be clear now that the freedom in 
redefining JF 7<5 — > T^s + fyS is equivalent to the changes r a — > r a + kn~ l b a we considered 
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earlier in this section. This freedom can be exploited together with the symmetries to 
simplify the problem of solving explicitly for r a in any specific MHD flow. 
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